Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandias DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to find optimal designs for an interdigitated flushwall injector suitable for scramjet applications at hypervelocity conditions. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. A Latin hypercube sampling design-of-experiments method was used to select design points for RAS. A methodology was developed that automated building geometries and generating grids for each design. The ensuing RAS analysis generated the performance database from which the two objective functions of interest were computed using a one-dimensional performance utility. The data were fitted using four surrogate models: an artificial neural network (ANN) model, a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model obtained an optimal solution set that predicted high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts are required in order to lower the errors and perform more accurate surrogate-based optimization. sed optimization

The Effect of Turbulence Modeling on the Mixing Characteristics of Several Fuel Injectors at Hypervelocity Flow Conditions

CFD analysis is presented on the effects of turbulence modeling choices on the mixing characteristics and performance of three fuel injectors at hypervelocity flow conditions. The analyses were carried out with the VULCAN-CFD solver using Reynolds-Averaged Simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The three injectors are the baseline configurations used in the experiments and represent three categories of injectors typically considered individually or in combination for fueling high-speed propulsive devices. The current work discusses the impact of the turbulence model and the turbulent Schmidt number on the mixing flow field behavior and the mixing performance as described by the one-dimensional values of the Mach number, total pressure recovery, and the mixing efficiency. Because planar laser induced fluorescence (PLIF) images are available from the EIMP experiments, the sensitivity of the synthetic LIF signal to turbulence modeling choices is also examined to determine whether PLIF can be extended beyond its intended qualitative visualization purpose and used to guide CFD turbulence model and parameter selections. It is found that the mixing performance, as quantified using mixing efficiency, exhibits a strong sensitivity to both turbulence model choice and turbulent Schmidt number value. However, the synthetic LIF signal only demonstrates a modest level of sensitivity, which suggests that PLIF is of limited use for guiding CFD turbulence model and parameter selections

Pre-Test CFD for the Design and Execution of the Enhanced Injection and Mixing Project at NASA Langley Research Center

With the increasing costs of physics experiments and simultaneous increase in availability and maturity of computational tools it is not surprising that computational fluid dynamics (CFD) is playing an increasingly important role, not only in post-test investigations, but also in the early stages of experimental planning. This paper describes a CFD-based effort executed in close collaboration between computational fluid dynamicists and experimentalists to develop a virtual experiment during the early planning stages of the Enhanced Injection and Mixing project at NASA Langley Research Center. This projects aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships relevant to flight Mach numbers greater than 8. The purpose of the virtual experiment was to provide flow field data to aid in the design of the experimental apparatus and the in-stream rake probes, to verify the nonintrusive measurements based on NO-PLIF, and to perform pre-test analysis of quantities obtainable from the experiment and CFD. The approach also allowed for the joint team to develop common data processing and analysis tools, and to test research ideas. The virtual experiment consisted of a series of Reynolds-averaged simulations (RAS). These simulations included the facility nozzle, the experimental apparatus with a baseline strut injector, and the test cabin. Pure helium and helium-air mixtures were used to determine the efficacy of different inert gases to model hydrogen injection. The results of the simulations were analyzed by computing mixing efficiency, total pressure recovery, and stream thrust potential. As the experimental effort progresses, the simulation results will be compared with the experimental data to calibrate the modeling constants present in the CFD and validate simulation fidelity. CFD will also be used to investigate different injector concepts, improve understanding of the flow structure and flow physics, and develop functional relationships. Both RAS and large eddy simulations (LES) are planned for post-test analysis of the experimental data

Comparison of Mixing Characteristics for Several Fuel Injectors at Mach 8, 12, and 15 Hypervelocity Flow Conditions

CFD analysis is presented of the mixing characteristics and performance of three fuel injectors at flight Mach numbers of 8, 12, and 15. The Reynolds-averaged simulations (RAS) were carried out using the VULCAN-CFD solver. The high Mach number flow conditions match those of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies relevant to flight Mach numbers greater than 8. The injectors include a fuel placement device, a strut, and a fluidic vortical mixer, a ramp. These fuel injectors accomplish the necessary task of distributing and mixing fuel into the supersonic cross-flow, albeit via different strategies. For comparison, a flush-wall injector is also included. This type of injector generally represents the simplest method of introducing fuel into a scramjet combustor. The three injectors represent the baseline configurations of the EIMP experiments. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are computed from the RAS and compared for the three flight conditions and injector configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations. Plotting the total pressure recovery and thrust potential as a function of mixing efficiency provides added insight into critical aspects of combustor performance as the flight condition and injector type are varied

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows

Comparison of Several Global Mixing Performance Metrics for High-Speed Fuel Injectors

To experimentally assess and compare the mixing performance of high-speed fuel injectors for scramjet engines, quantitative global metrics are needed. The one-dimensional metric most commonly used to assess the degree of mixing completeness at a given downstream station is the mixing efficiency parameter. The experimental determination of the mixing efficiency parameter requires measurement of the spatial distributions of both the fuel mass fraction and the mass flux. Standard in-stream gas sampling techniques can be used to measure the fuel mass fraction distribution, however the mass flux distribution is not easily determined experimentally because it requires the measurement of three independent aerothermodynamic variables in addition to the mixture composition. For this reason, several metrics that can be calculated from the fuel distribution alone are commonly used to assess mixing performance. Because these other metrics do not provide a mass flux-weighted measure of the local degree of mixing completeness, they may not correlate well with the mixing efficiency parameter. Therefore, if the substitute metrics are to be used to compare the mixing performance of candidate fuel injector concepts, it is important to understand their relationships to the mixing efficiency parameter in a representative scramjet combustor flowfield. This work investigates the relationships between the mixing efficiency parameter and several substitute metrics that are able to be measured with the current experimental setup of the Enhanced Injection and Mixing Project at the NASA Langley Research Center for baseline strut and ramp injectors. The results of these comparisons have revealed that it is possible to glean different (i.e., incorrect) conclusions about which injector is the better mixer when the substitute mixing performance metrics are used instead of the mixing efficiency parameter, thereby highlighting the importance of mass flux-weighted mixing performance metrics

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach

A Priori Analysis of a Compressible Flamelet Model using RANS Data for a Dual-Mode Scramjet Combustor

In an effort to make large eddy simulation of hydrocarbon-fueled scramjet combustors more computationally accessible using realistic chemical reaction mechanisms, a compressible flamelet/progress variable (FPV) model was proposed that extends current FPV model formulations to high-speed, compressible flows. Development of this model relied on observations garnered from an a priori analysis of the Reynolds-Averaged Navier-Stokes (RANS) data obtained for the Hypersonic International Flight Research and Experimentation (HI-FiRE) dual-mode scramjet combustor. The RANS data were obtained using a reduced chemical mechanism for the combustion of a JP-7 surrogate and were validated using avail- able experimental data. These RANS data were then post-processed to obtain, in an a priori fashion, the scalar fields corresponding to an FPV-based modeling approach. In the current work, in addition to the proposed compressible flamelet model, a standard incompressible FPV model was also considered. Several candidate progress variables were investigated for their ability to recover static temperature and major and minor product species. The effects of pressure and temperature on the tabulated progress variable source term were characterized, and model coupling terms embedded in the Reynolds- averaged Navier-Stokes equations were studied. Finally, results for the novel compressible flamelet/progress variable model were presented to demonstrate the improvement attained by modeling the effects of pressure and flamelet boundary conditions on the combustion

Comparison of Mixing Characteristics for Several Fuel Injectors on an Open Plate and in a Ducted Flowpath Configuration at Hypervelocity Flow Conditions

In order to reduce the cost and complexity associated with fuel injection and mixing experiments for high-speed flows, and to further enable optical access to the test section for nonintrusive diagnostics, the Enhanced Injection and Mixing Project (EIMP) utilizes an open flat plate configuration to characterize inert mixing properties of various fuel injectors for hypervelocity applications. The experiments also utilize reduced total temperature conditions to alleviate the need for hardware cooling. The use of "cold" flows and non-reacting mixtures for mixing experiments is not new, and has been extensively utilized as a screening technique for scramjet fuel injectors. The impact of reduced facility-air total temperature, and the use of inert fuel simulants, such as helium, on the mixing character of the flow has been assessed in previous numerical studies by the authors. Mixing performance was characterized for three different injectors: a strut, a ramp, and a flushwall. The present study focuses on the impact of using an open plate to approximate mixing in the duct. Toward this end, Reynolds-averaged simulations (RAS) were performed for the three fuel injectors in an open plate configuration and in a duct. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are then computed and compared for the two configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations